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^ Abstract 

b ' 

One of the outstanding problems of the dynamical evolution of the outer solar system concerns the observed population ratio 
between the Oort Cloud (OC) and the Scattered Disc (SD): observations suggest that this ratio lies between 100 and 1 000 but 
simulations that produce these two reservoirs simultaneously consistently yield a value of the order of 10. Here we stress that the 
populations in the OC and SD are inferred from the observed fluxes of new Long Period Comets (LPCs) and Jupiter-family comets 
(JFCs), brighter than some reference total magnitude. However, the population ratio estimated in the simulations of formation of 
i - 1 the SD and OC refers to objects bigger than a given size. There are multiple indications that LPCs are intrinsically brighter than 
JFCs, i.e. an LPC is smaller than a JFC with the same total absolute magnitude. When taking this into account we revise the 
l~M SD/JFC population ratio from our simulations relative to Duncan and Levison (1997), and then deduce from the observations that 
r"| the size-limited population ratio between the OC and the SD is 44+^. This is roughly a factor of four higher than the value 12 + 
Qh!1 that we obtain in simulations where the OC and the SD form simultaneously while the planets evolve according to the so-called 
q 'Nice model' . Thus, we still have a discrepancy between model and 'observations, but the agreement cannot be rejected by the null 
hypothesis. 
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^ 1. Introduction and background 

OC ■ 

When examining the orbital data of new comets entering 
the inner solar system, Oort (1950) discovered that the distri- 
ct bution of reciprocal semi-major axis of these comets showed a 
£C) distinct excess for I /a < 5 x 10~ 4 ALT 1 . The observed semi- 
C> major axis distribution led Oort to suggest that the Sun is sur- 
rounded by a cloud of comets in the region between 20 000 AU 
I to 150000 AU, and that it contains approximately 10 11 comets 
J> with isotropic inclination and random perihelia. This hypoth- 
k> esised cloud of comets surrounding the Sun is now called the 
%J "Oort cloud' (OC). 

.9: 

The formation and evolution of this reservoir of comets has 
been an issue of study ever since its discovery. The main uncer- 
tainties are its population, how it formed and how its existence 
ties in with what we know about the evolution of the rest of the 
outer solar system. We review each of these below and then 
state the motive behind this study. 
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1. 1. Total population of the Oort cloud 

The only method with which we can infer the number of 
comets in the OC is by determining the flux of long-period 
comets (LPCs), which can be divided into new comets (NCs) 
and returning comets (RCs). The NCs are a proxy for the total 
population of the OC: the total population of the cloud can be 
inferred from their flux through the inner solar system (Wiegert 
& Tremaine, 1999). New comets are traditionally classified 
as those with semi-major axis a > 10 kAU (e.g. Wiegert & 
Tremaine, 1999). 

There exist two agents which perturb the comets in the cloud 
onto orbits that enter the inner solar system: passing stars (Weiss- 
man, 1980; Hills, 1981) and the Galactic tide (Heisler & Tremaine, 
1986; Levison et al., 2001). The passing stars cause usually 
small random deviations in the orbital energy and other orbital 
elements. The Galactic tide on the other hand systematically 
modifies the angular momentum of the comets at constant or- 
bital energy. Heisler & Tremaine (1986) discovered that if the 
semi-major axis of the comet is large enough then the comet's 
change in perihelion can exceed 10 AU in a single orbit and 
'jump' across the orbits of Jupiter and Saturn and thus not suf- 
fer their perturbations before becoming visible; in this partic- 
ular case the comet is considered an NC. However, Heisler 
& Tremaine's (1986) approximation of the Galactic tide only 
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used the vertical component, which is an order of magnitude 
stronger than the radial components. In this approximation the 
z-component of the comet's orbital angular momentum in the 
Galactic plane is conserved and the comets follow closed tra- 
jectories in the q - to plane (with q being the perihelion dis- 
tance and u> being the argument of perihelion). Including the 
radial tides breaks this conservation and the flux of comets to 
the inner solar system from the OC is increased (Levison et al., 
2006). The trajectories that lead comets into the inner solar sys- 
tem should be quickly depleted, were it not for the passing stars 
to refill them (Rickman et al, 2008). The synergy between these 
two perturbing agents ensures there is a roughly steady supply 
of comets entering the inner solar system. 

Even though the perturbations on the cloud are now under- 
stood, there remain large uncertainties in the total number of 
comets in the cloud and there have been many attempts to con- 
strain it (e.g. Oort, 1950; Hills, 1981; Weissman, 1983, 1996; 
Weigert & Tremaine, 1999; Francis, 2005). The general con- 
sensus seems the total number is between 10 11 to 10 12 comets 
with total absolute magnitude Hj < 1 1 . Here Hj is given by 

H T = V dls -51og£>-2.5vlogr (1) 

where Vdi s is the apparent magnitude (nucleus with coma), D 
is the distance of the comet to Earth, v is a measure of how the 
brightness scales with heliocentric distance (the photometric in- 
dex), and r is the distance of the comet to the Sun. Thus, Hj is 
the magnitude of a comet (nucleus plus coma) if viewed from 
the Sun, at a distance of 1 AU. There is a great variability in v 
among comets. The first to catalogue this quantity for a large 
sample of comets was Whipple (1978), who found that for NCs 
on their inbound leg (v) = 2.44 + 0.3 and on their outbound leg 
(v) = 3.35 ± 0.4. For short-period comets v is usually higher 
than 3. From a limited sample of LPCs Sosa & Fernandez 
(2011) find (v) ~ 3. Using v = 4 yields to the commonly- 
used value Z/io, which is close to Whipple's (1978) average for 
both LPCs and short-period comets. A value Hj < 1 1 is used 
throughout the literature when referring to the total number of 
comets in the OC, and thus we do so as well. 

The best estimate of the number of comets in the cloud 
comes from a dynamical study by Kaib & Quinn (2009). Based 
on a suggestion by Levison et al. (2001), Kaib & Quinn (2009) 
find that the fraction per unit time of the OC that is visible in 
the form of a new comet is 10" 11 yr _1 . We shall adopt this 
fraction at a later stage in this paper. By calibrating this prob- 
ability to the flux of new comets from Francis (2005) Kaib & 
Quinn (2009) conclude that the whole OC contains (2-3)xlO n 
comets. Thus the most recent simulations, and observations, 
suggest that the cloud contains approximately (1 - 5) x 10 n 
comets, rather than 10 12 , of size equivalent to that of an LPCs 
with Hj < 1 1 . 

1.2. Formation and evolution 

The first attempt to form the OC by direct numerical inte- 
gration was undertaken by Duncan et al. (1987), who found that 



comets with original perihelion q > 15 AU were likely to reach 
the cloud while those with shorter perihelion distance were not. 
Duncan et al. (1987) found that the inner edge of the cloud is 
located at approximately 3 000 AU while the (assumed) outer 
edge is at 200 000 AU. However, some aspects of their results, 
in particular their high formation efficiency, are questionable 
because their initial conditions turn out not to be representative 
of the reality: their assumption that the first stage of scattering 
does not greatly change the perihelion distance is incorrect. 

Dones et al. (2004) performed a study similar to Duncan et 
al. (1987) but with more realistic initial conditions. At the end 
of the simulation Dones et al. (2004) obtained a formation effi- 
ciency of only 5%. Similar results are reported elsewhere (e.g. 
Kaib & Quinn, 2008; Brasser et al., 2010). 

However, simulations of OC formation have suffered from a 
difficult problem: they are unable to reproduce the inferred ratio 
between the population of the Scattered Disc (SD) (Duncan & 
Levison, 1997) and the OC. The SD is believed to be the source 
of the Jupiter-family comets (JFCs) (Duncan & Levison, 1997). 
The JFCs are a set of comets whose Tisserand parameter with 
respect to Jupiter satisfies Tj e (2,3] (Levison, 1996). Care 
has to be taken here because this criterion would erroneously 
classify some comets as JFCs when they have Tj < 3 due to i) 
a high inclination but are entirely outside of Jupiter's orbit, or 
ii) very low inclination objects that have their perihelia close to 
Jupiter but very long semi-major axis (Gladman et al., 2008). 
Thus we follow Gladman et al. (2008) and additionally impose 
that a JFC must have q < 7.35 AU so that the object is most 
likely dynamically controlled by Jupiter. 

From the current population of JFCs, Duncan & Levison 
(1997) and Levison et al. (2008a) estimate there are 6 x 10 8 
bodies in the SD whose size is equivalent to that of JFCs with 
Hj < 9. This estimate of the SD population is approximately a 
factor of 150 to 750 lower than the estimate of the OC popula- 
tion. Compensating for the difference in Hj used in the deriva- 
tions for the OC (H T < 11) and SD (H T < 9) populations is 
difficult because JFCs and LPCs appear to have strongly differ- 
ent //-distributions in the 9-11 range (Fernandez et al., 1999; 
Fernandez & Sosa, 2012; Francis, 2005). Numerical simula- 
tions tend to yield a 5:1 (Kaib & Quinn, 2008) to 20:1 ratio 
(Dones et al, 2004; Brasser et al., 2010). Thus, simulations 
consistently underestimate the observed OC to SD population 
ratio by at least an order of magnitude! This large discrepancy 
forms the motivation of our study. 

Two solutions have been proposed to remedy this problem: 
forming part of the OC while the Sun was still in its birth cluster 
(Fernandez & Brumni, 2000; Brasser et al., 2006, 2012; Kaib 
& Quinn, 2008), or forming the OC by the capture of comets 
from other stars (Levison et al., 2010). Even though these sce- 
naria appear to be able to solve the SD to OC population dis- 
crepancy, unfortunately both could suffer from the difficulty of 
scattering small comets to large distances in the presence of gas 
drag (Brasser et al., 2007). In light of the above problems, we 
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decided to re-examine the whole problem from scratch, as de- 
tailed below. 

1.3. Our approach 

Scattering small comets to large distances in the presence of 
gas is very difficult. Therefore we suggest that the SD and the 
OC formed together after the removal of the gas in the proto- 
planetary disc. The Nice model provides the natural framework 
for such a 'late' (relative to gas removal) and contemporary for- 
mation of both these reservoirs. The Nice model argues that, 
during the gas-disc phase, the giant planets had orbits more cir- 
cular and were more closely packed than now. The current plan- 
etary orbits have been achieved during a phase of dynamical in- 
stability of the planets that occurred after gas removal (Tsiganis 
et al., 2005; Morbidelli et al., 2007), possibly as late as the Late 
Heavy Bombardment event (approximately 500 Myr after gas 
removal; Gomes et al., 2005; Levison et al., 201 1). During this 
planetary instability, a primordial trans-Neptunian disc of plan- 
etesimals was dispersed, with just a few of its objects surviving 
today in the Kuiper Belt, in the SD, in the OC and in the Tro- 
jan populations of Jupiter and Neptune (Morbidelli et al., 2005; 
Levison et al. 2008a,b). The previous simulations addressing 
the OC/SD ratio (e.g. Dones et al. 2004,; Kaib & Quinn, 2008; 
Brasser et al., 2010) assumed that the planets were on their cur- 
rent orbits, which is an unlikely scenario. Thus, in section 4. 1 
we re-asses this ratio in the framework of the Nice model, where 
the different orbital evolution of the planets might in principle 
lead to a different result. 

The higher eccentricities and inclinations of Uranus and 
Neptune during their migration produces a dynamically hot SD, 
where the inclination distribution ranges up to several tens of 
degrees, consistent with the observed distribution of large SD 
objects. This SD has a very different structure from that used in 
Duncan & Levison (1997), which was dynamically cold (incli- 
nations up to approximately 15°). It is widely believed that the 
JFCs originate in the SD and thus the two should be intimately 
linked (Duncan & Levison, 1997). In principle, the population 
ratio between the SD and JFCs after giant planet migration may 
be different from that estimated in Duncan & Levison (1997). 
We investigate this in section 4.2, leading to a new estimate 
of the SD population from the observed JFC population. The 
OC produced in the Nice model is equivalent to that formed in 
previous studies (e.g. Dones et al., 2004; Kaib & Quinn, 2008; 
Brasser et al., 2010), and therefore we do not need to re-evaluate 
the relationship between the OC population and the LPC flux. 

Finally, in section 5 we re-examine the LPC and JFC popu- 
lations. We take into account that there are multiple indications 
that LPCs are intrinsically brighter than JFCs: an LPC of size 
comparable to that of a JFC with Hj ~ 9 has Hj ~ 6.5. To- 
gether with the results of section 4.2 we re-assess the 'observed' 
OC/SD ratio for a size-limited population. These are then com- 
pared with the results obtained from the Nice-model simulation 
from section 4. 1 . Finally, we draw our conclusions in section 6. 
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Figure 1 : Evolution of the planets. The curves correspond to semi-major axis, 
pericentre and apocentre of the planets. Jupiter stays close to 5 AU, Saturn 
close to 9 AU, Uranus ends up close to 19 AU and Neptune settles at 3 1 AU. 



2. Planetary evolution 

For the purpose of this study we have used one of the Nice- 
model planetary evolutionary tracks presented in Levison et al. 
(2008b). More precisely, using the recipe described in Levison 
et al. (2008b), we re-enact the evolution shown in their Run 
A in which, at the end of the phase of mutual close encounters 
among the planets, Neptunes semi-major axis is a^ = 27.5 AU 
and its eccentricity is = 0.3. Uranuss semi-major axis and 
eccentricity are ajj = 17.5 AU and en = 0.2. The mutual in- 
clination of both planets is approximately 1 degree. The evo- 
lution of the planets after the mutual encounters is depicted in 
Fig. [T] where we plot the semi-major axis, perihelion distance 
and aphelion (Q) distance for the four giant planets as a func- 
tion of time. Jupiter remains at 5 AU, Saturn stays at 9.5 AU. 
Uranus migrates from 17 AU to nearly 19 AU while Neptune 
migrates outwards to settle close to 31 AU. 

We want to emphasize that the real evolution of the planets 
cannot be traced, so that we do not expect that the evolution 
we consider is exactly right. However, the evolution of Run 
A leads to final planetary orbits very similar to the current ones 
and shows a high compatibility with the currently known orbital 
structure of the Kuiper Belt (Levison et al., 2008b). Hence we 
argue the evolution above is representative of what could have 
happened in reality. 

Having decided on the evolution of the planets we explain 
our numerical methods in the next section. 



3. Numerical methods 

In this section we describe our numerical methods. We ran 
two series of simulations: one set for the formation of the OC 
and another set for the formation and evolution of the SD. The 
planetary evolution for both sets of simulations is the same. The 
reason we examine the SD separately is because we need a high 



3 



resolution (large number of comets) to determine the JFC to 
SDO population ratio. This high number of comets in the SD 
simulations was achieved by repeated cloning of the remaining 
comets at several stages. In addition, for the evolution of the 
SD, in particular for SDOs to become JFCs, we do not need the 
influence of the Galactic tide and passing stars. 

3.1. Oort cloud 

The simulations that were performed for the formation of 
the OC consist of two stages. During the first stage, the plane- 
tary evolution shown in Fig. Q] is re-enacted for a total duration 
of 200 Myr. We added 6 000 massless test particles to each sim- 
ulation. Their initial conditions were taken from Levison et al. 
(2008b): the semi-major axes were between 29 AU and 34 AU, 
their eccentricities were 0.15 and their orbits were coplanar. 
The time-step for the simulations was 0.4 yr. The perturba- 
tions from the Galactic tide were included using the method of 
Levison et al. (2001), which incorporates both the vertical and 
radial Galactic tides. The local Galactic density was 0.1 M Q 
pc -3 (Holmberg & Flynn, 2000). We assumed a flat Galactic 
rotation curve with the Sun having an angular velocity of 30.5 
km s kpc 1 (MacMillan & Binney, 2010). Passing stars were 
included according to the method described in Rickman et al. 
(2008) at a Galactic distance of 8 kpc, with the Sun's sphere of 
influence being 1 pc. We ran a total of five realisations. The 
difference between each simulation lies in the different initial 
conditions for the test particles. 

After the first stage was completed, we took the positions 
and velocities of the planets and test particles and resumed the 
integration for another 4 Gyr using SCATR (Kaib et al., 201 1), 
which is a Symplectically-Corrected Adaptive Timestepping Rou- 
tine. It is based on SWIFT's RMVS3 (Levison & Duncan, 
1994). It has a speed advantage over SWIFT's RMVS3 or 
MERCURY (Chambers, 1999) for objects far away from both 
the Sun and the planets where the time step is increased. We set 
the boundary between the regions with short and long time step 
at 300 AU from the Sun (Kaib et al. 201 1). Closer than this dis- 
tance the computations are performed in the heliocentric frame, 
like SWIFT's RMVS3, with a time step of 0.4 yr. Farther than 
300 AU, the calculations are performed in the barycentric frame 
and we increased the time step to 50 yr. The error in the energy 
and angular momentum that is incurred every time an object 
crosses the boundary at 300 AU is significantly reduced through 
the use of symplectic correctors (Wisdom et al., 1996). For the 
parameters we consider, the cumulative error in energy and an- 
gular momentum incurred over the age of the solar system is of 
the same order or smaller than that of SWIFT's RMVS3. The 
same Galactic and stellar parameters as in the first simulation 
were used. Comets were removed once they were further than 
1 pc from the Sun, or collided with the Sun or a planet. 

3.2. Scattered Disc 

For the simulations of the SD we used the following strat- 
egy. We took the planets and comets at the end of the planetary 



migration phase and removed all comets that were further than 
3 000 AU from the Sun. Unlike the above case for the OC, to 
correctly simulate the decay of the SD the planets need to be 
on their current orbits, or match these as closely as possible. 
Therefore after the first stage of migration was completed, we 
performed a second stage where we artificially migrated Uranus 
outwards by 0.25 AU to its current orbit over a time scale of 

5 Myr. We kept Neptune's semi-major axis at its final position 
of 30.7 AU because lowering it to its current value of 30. 1 AU 
would cause resonant objects in the SD to escape from their res- 
onances. We also had to damp the eccentricities of Uranus and 
Neptune somewhat. Neptune's eccentricity was lowered from 
0.015 to 0.01 while Uranus' eccentricity was lowered from 0.06 
to its current value of 0.04. However, it matches its current 
secular properties and at its secular maximum it is 0.06. The 
migration was accomplished by interpolation of the planets' or- 
bital elements with SWIFT RMVS3 (Petit et al., 2001; Brasser 
et al., 2009; Morbidelli et al., 2010). During this migration we 
kept the precession frequencies of the planets as close as pos- 
sible to their current values. At the end of this fictitious migra- 
tion and eccentricity damping the giant planets had their current 
secular architecture. During this short simulation comets were 
removed once they hit a planet, hit the Sun or were farther than 
3 000 AU from the Sun. The passing stars and the Galactic tide 
were not included. 

After Uranus was artificially migrated and the planets resided 
on their correct orbits, we proceeded to integrate the planets 
and comets for another 3.8 Gyr using SWIFT RMVS3. We 
cloned all comets that remained after the migration three times. 
Cloning was achieved by adding a random deviation of 10 -6 
radians to the comets' mean anomaly, keeping all the other ele- 
ments fixed. We stopped the simulations at 1 Gyr and 3.5 Gyr to 
clone the remaining comets three times. This repeated cloning 
ensured enough comets remained during the last 500 Myr for 
good statistics on JFC production. 

During the last 500 Myr it was essential to keep track of 
visible JFC production since we shall use this population as a 
proxy for the number of comets in the Scattered Disc (Duncan 

6 Levison, 1997). Here we copy the term 'visible JFC from 
Levison & Duncan (1997) to refer to a JFC with perihelion dis- 
tance q < 2.5 AU. We modified SWIFT RMVS3 to output all 
particles that have q < 2.5 AU every 100 yr. Afterwards we fil- 
ter out the visible JFCs by requiring they obey q < 2.5 AU 
and Tj e [2,3]. Here Tj is the Tisserand parameter of the 
comet with respect to Jupiter. We also ran the last 500 Myr 
in a second set of simulations where we emulated Levison & 
Duncan (1997) and removed any particle the moment it came 
closer than 2.5 AU from the Sun. These last simulations were 
done to accurately determine the fraction of SDOs that become 
visible JFCS, / v jfc- 

All simulations were performed on either TIARA Grid or 
the ASIAA Condor pool. For the OC each simulation lasted 
only about a week, thanks to SCATR's speed. For the SD simu- 
lations, reaching 4 Gyr took a couple of months of computation 
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Figure 2: Top panel: Percentage of comets trapped in the OC. The bottom panel 
depicts the population ratio between the OC and the SD. 



time per simulation. 

4. Results from numerical simulations 

4.1. Formation of the Oort Cloud and Scattered Disc and pre- 
diction of the OC/SD ratio 
There have been many several publications that have stud- 
ied OC formation in the current galactic environment (Dones 
et al., 2004; Kaib & Quinn, 2008; Dybczyriski et al., 2008; 
Neslusan et al., 2009; Leto et al., 2009; Brasser et al., 2010) 
and thus we choose not to do an in-depth analysis of the struc- 
ture of the cloud. Instead we list a few key issues because most 
properties of the cloud produced in the Nice model are unlikely 
to be very different from previous works. 

We define a comet to be in the OC when it has both a > 
1 000 AU and q > 40 AU because when a ~ 2000 AU the 
Galactic tide begins to dominate over planetary perturbations 
from Neptune. The condition q > 40 AU is imposed to en- 
sure that the object has been decoupled from Neptune by the 
Galactic tide. A SD object is defined to have a < 1 000 AU and 
q > 30 AU because its motion is controlled by Neptune. For the 
sake of completeness, an object with q e [5, 30] AU and not in 
resonance with Neptune is considered a Centaur and an object 
with q e [30,40] and a > 1 000 AU is called a high-a SDO 
(HaSDO) even though these could be low-q OC objects. For 
all three populations no restrictions were placed on the inclina- 
tion. We realise that our classification is not complete because 
it leaves out resonant Neptune-crossing objects like Pluto, but 
since these only comprise a very small subset of all objects we 
believe our classification is justified. Thus, we restrict ourselves 
to a < 1 000 AU for SD objects (SDOs). 

In Fig. |2] we have plotted the percentage of the original 
comets that are part of the OC (top panel) and the population 
ratio between the OC and the SD (bottom panel) for one sim- 
ulation. The results from the other simulations are very simi- 
lar. We find that the average OC formation efficiency at 4 Gyr 



is 7.1% + 0.3% (1-sigma error, used throughout this paper). 
This value is higher than that reported elsewhere in the litera- 
ture (e.g. Dones et al. 2004; Kaib & Quinn 2008; Dybczyhski 
et al. 2008; Brasser et al., 2010). The reason for the higher effi- 
ciency is twofold. First, all of our comets were initially placed 
in the vicinity of Neptune. Second, Uranus and Neptune were 
initially on more eccentric orbits so that they are capable of 
keeping more objects under their dynamical control rather than 
passing them down to Saturn and Jupiter. It is well-known that 
Uranus and Neptune place objects into the OC while Jupiter and 
Saturn tend to eject them (e.g. Duncan et al., 1987). 

From our simulations we deduce that the population ratio 
between the OC and the SD is 12 + 1. This ratio is a little lower 
than that found by Dones etal. (2004) and Brasser etal. (2010), 
but similar to that of Kaib & Quinn (2008). We attribute this to 
our higher OC formation efficiency. 

Having summarised the key properties of the OC we now 
turn to the SD and JFC production. 

4.2. Linking the SD and JFC populations 

Before addressing the population ratio between the SD and 
the JFC population we need to demonstrate that the SD pro- 
duced in our simulations is an acceptable source for the JFCs. 
In other words, we need to show that the the objects that evolve 
from the SD into the JFC region reproduce the observed distri- 
bution of the actual JFCs acceptably well. 

Here we use the method described in Levison & Duncan 
(1997), who followed a large number of test particles from 
the Kuiper Belt into the visible JFC region. Levison & Dun- 
can (1997) show that the inclination distribution of their simu- 
lated visible JFCs on their first apparition is inconsistent with 
the observed one. Knowing that the typical inclination of the 
JFCs increases with time they use the real distribution to deter- 
mine the physical lifetime of the comets (also referred to as the 
fading time). They compare the inclination distributions from 
their simulations with the observed one for various values of 
the physical age of the comets. 

To do this, Levison & Duncan (1997) define £(r,i) as the 
number of comets with inclinations between i and i+di and with 
physical ages between t and r+dr. Here t is measured from the 
time when the comet first enters into the visible (q < 2.5 AU) 
region. By assuming that all comets fade instantaneously after 
a certain time Tf and remain dormant afterwards, the inclination 
distribution of active comets is thus given by 



6(0 



/ 



£(t, i)dr. 



(2) 



We define similar distributions for the Tisserand parameter, £; d (Tj), 
and the minimum distance between Jupiter's orbit and one of 
the comet's nodes, ^ a (dj). Here dj = min(|ay - ra\, \dj - ry|), 
where rn.u - a (l _ g2 )/(l + e cos at) are the distances of the 
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Figure 3: Top-left panel: Probability that the inclination (red), dj (blue) and Tj 
(black) distributions of the visible JFCs produced in our simulations matches 
their observed ones as a function of the ageing time, t. Top-right panel: the 
median value of dj from our simulations as a function of ageing time. The 
red line indicates the observed value. Bottom-left: the ratio of the number of 
dormant to active 'visible' JFCs as a function of ageing time. Bottom-right: the 
median value of Tj. The red line indicates the observed value. 

comet's nodes to the Sun. Most JFCs are believed to be scat- 
tered onto low-g orbits by encounters with Jupiter and thus the 
distribution of dj should be indicative of the dynamical age of 
these comets (Levison & Duncan, 1997). Levison & Duncan 
(1997) obtain £,(/) and /) for a range of Tf from their simu- 
lations. They compute the probability of a match between their 
simulations and the real comets by performing a Komolgorov- 
Smirnov (KS) test for various values of Tf. The peak of the KS 
probability as a function of the active lifetime t yields the best- 
fit value of Tf . Their distribution peaks at Tf = 12 kyr and they 
conclude that this value must correspond to the typical fading 
time for the visible JFCs. 

We have repeated their procedure but included also the re- 
sults from ^(Tj). The outcome is shown in Fig. [3] In the top- 
right panel we depict the probability of the inclination (red), d j 
(blue) and Tj (black) distributions matching their observed dis- 
tributions as a function of the physical lifetime, t. From the plot 
it appears that the best-fit value for the physical lifetime, Tf , lies 
between 10 kyr and 15 kyr, bracketing the value of 12 kyr found 
by Levison & Duncan (1997). Thus, we also adopt Tf = 12 kyr 
here. Our match between the inclination and d j distribution is 
not as good as Levison & Duncan (1997) nor is the peak of the 
latter as high, but it is nevertheless a satisfactory match (KS 
probability larger than 20%). The top-right panel depicts the 
median dj as a function of t, with the red line indicating the 
observed value. The corresponding distribution for Tj is shown 
in the bottom-right panel. Note that the median d j and Tj dis- 
tributions cross their observed value at different t, indicating 
that the dynamics is not entirely controlled by scattering off of 
Jupiter. Last, the bottom-left panel of Fig. [3] shows the ratio 
of the total number of JFCs (active and dormant) to the num- 
ber of active ones. The total number of comets with t < Tf is 
computed as r/(Tf) = J £; R (i)di and the total number of JFCs is 



Figure 4: Plot of \ ja vs Tj of the observed (open circles) and simulated comets 
(bullets) with q < 2.5 AU and with Tj < 3.1. The data from the simulations are 
for the duration of the comets' active lifetime. 

computed by taking t — > oo. 

For a nominal Tf = 12 kyr, the total number of comets is 
about 6.5 times the number of active ones, similar to the factor 
5 what was reported in Levison & Duncan (1997) and Di Sisto 
et al. (2009). 

There is an additional check we performed. We checked 
whether we generated a Saturn-family comets group i.e. comets 
that were scattered into the inner Solar System by Saturn and 
have a node and/or their aphelion near Saturn. This would be 
problematic, given that these comets are not observed. The 
Saturn-family comets would have Tj < 2.6 and a > 6 AU. In 
Fig.|4]we have plotted Tj vs 1 /a for all comets with P < 200 yr, 
q < 2.5 AU and Tj < 3.1 (open circles). The observational 
data for this plot was obtained from jpiB We superimposed 
the comets from our simulations during their active lifetime 
(bullets). This lifetime is usually 12 kyr but shorter for some 
comets who were dynamically removed from the visible region 
before having faded. As one may see the distribution of bullets 
with (I /a, Tj) — (0.15,2.6) and lower are compatible with the 
real comets. Thus we do not generate a swarm of Saturn-family 
comets. 

We performed a few further tests, such as establishing the 
fraction of visible comets that are Halley-types, the average 
time spent as a Centaur, the visible JFC perihelion distribution 
and the average dynamical lifetime of JFCs. The former two 
agree well with the results of Levison & Duncan (1997), but we 
find a shorter mean JFC dynamical lifetime: 165 kyr + 60 kyr 
vs. 270 kyr for Levison & Duncan (1997). 

In summary, we have a working model linking the SD pop- 
ulation to the visible JFCs population and a reliable estimate for 
the physical lifetime of comets. We can now use this informa- 



'http://ssd.jpl.nasa.gov/sbdb_query.cgi 
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Figure 5: Remaining fraction of comets in the SD. The result of five simulations 
are shown. The solid line shows the remaining fraction from from Duncan & 
Levison (1997). Data kindly provided by Hal Levison. The dashed lines shows 
the results our simulations. 



tion to compute the SDO to JFC population ratio. 

The formula relating the number of objects in the SD (A^do) 
to the number of active objects in the JFC visible region (A^jfc) 
is: 

AW 



A^SDO - 



1"vJFcksDl/vJFC 



(3) 



where tsd is the fractional decay rate of the SD at the current 
time, /vJfc is the fraction of the comets escaping from the SD 
that penetrate into the visible region and t v jfc is the mean life- 
time spent by these comets in the visible JFC region as active 
comets. We evaluate these three quantities below. 

4.2.1. Evaluation of the fractional decay rate of the SD, tsd 

An SDO is considered to have left the SD if it has been 
removed from the simulation (through ejection or collision) or 
if it is still evolving but has spent any time as a Centaur (i.e. 
has achieved q < 30AU). Defining by /sd(0 the fraction of 
the original trans-Neptunian disk surviving in the SD at time 
t (see Fig. 0, the fractional decay rate of the SD population 
is defined as rso(t) = (c//sd(0/^0//sd(0- We measured tsd 
over the last 0.5 Gyr of our simulations and obtained (/"sd) = 
-(1 .63 + 0.66) x 10~ 10 per year. For reference, Duncan & Lev- 
ison (1997) report a value of (rso) = _ 2.7 x 10~ 10 per year. 
This roughly factor of two difference is due to a higher fraction 
of our SD residing in a fossilized state (resonant and detached 
objects) than Duncan & Levison (1997). In fact, Gomes et al. 
(2005b) and Gomes (2011) have shown that the migration of 
Neptune can create a Fossilised Scattered Disc (FSD), which 
is comprised of planetesimals that are either no longer inter- 
acting with Neptune, such as (136199) Eris, 2000 CR105 and 
2004 XR190 (Buffy), or planetesimals that are trapped in mean- 
motion resonances with Neptune for a long time. 



4.2.2. Evaluation of the fraction ofJFCs, / v jfc 

To evaluate the fraction of SDOs that become visible JFCs 
precisely, we need to make sure we do not miss any object that 
enters the q < 2.5 AU region, even for a very short time. For this 
reason we have done an extra set of simulations where we re- 
move objects when they achieve this perihelion threshold. The 
check on the perihelion distance is done at every simulation 
timestep. 

We find that of the comets that leave the SD, an average 
fraction / v jfc = 16.5% + 8.0% penetrate into the region of vis- 
ible JFCs. Our value of / v jfc is somewhat lower than that re- 
ported in Levison and Duncan (1997) (30%). We attribute the 
difference to our source population having a different orbital 
distribution than that of Levison & Duncan (1997). 

4.2.3. Evaluation of the active visible lifetime, t v jfc 

The average time any JFC spends with q < 2.5 AU as an 
active object is computed using the simulations with a 100 yr 
high-resolution output. More precisely, we assumed that at each 
output entry the comet spends the entire 100 yr output time in 
the visible region. The comet is discarded once it reaches a 
physical lifetime of 12 kyr. We obtain t v jfc = 2.6 kyr, some- 
what shorter than the 3.6 kyr from Di Sisto et al. (2009). 

4.2.4. The Scattered Disc to visible JFC population ratio 

To compute the total number of SDOs we assume that the 
distributions in tsd and / v jfc are Gaussian with means and stan- 
dard deviations given by the nominal and error values listed 
above. Plugging the values of all of the variables and their 
distributions reported above into equation d3) we find a mean 
value of A^sdo = l-5*o'g X 10 7 A\,jfc- The error-bars are once 
again 1-sigmaj values. By comparison, this relationship yields 
N SD0 = 6.0 x 10 6 JV v jfc for Duncan & Levison (1997). We 
use this new SDO to JFC ratio in the next section to infer the 
OC/SD ratio from observations of size-limited populations of 
LPCs and JFCs. 

5. A re-examination of the Oort cloud to Scattered Disc 
population ratio 

In this section we attempt to obtain an updated value of the 
OC to SD population ratio for size-limited populations, rather 
than for populations limited in total absolute magnitude. In 
other words we want to obtain a value for the number of Oort 
cloud objects, Noc, and the number of SDOs, A^sd, for comets 
of the same size. We use the flux of LPCs in the inner solar 
system as a proxy for Noc (Wiegert & Tremaine, 1999) and 
the number of visible JFCs as a proxy for the number of SDOs 
(Duncan & Levison, 1997). The overall result depends on the 



2 The 1-sigma error values were computed with Monte Carlo simulations. 
We generated 10 5 values of / v jfc and »"sd assumed they followed a Gaussian 
distribution with the mean and the standard deviation reported above. From 
these we evaluated Nsu- The 1-sigma uncertainty interval for N^o is given by 
the values that leave 15.8% of the cumulative distribution of JVsd in each of the 
wings. 
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flux of LPCs that enter the inner solar system and on the ab- 
solute magnitudes of an LPC and a JFC of the same size. We 
discuss each of these below. 

5.1. The flux of new comets entering the inner solar system 

The flux of new comets entering the inner solar system is 
poorly known. The most-cited sources estimating the flux are 
Everhart (1967), Hughes (2001) and Francis (2005). Most of 
these list the flux of new comets with a perihelion distance 
q < 4 AU and absolute magnitude Hj < 10.9. However, there 
is a potential problem with using the flux of new comets with 
Hj < 10.9: the sample may be incomplete for new comets with 
Hj > 6.5. All three sources (Everhart, 1967; Hughes, 2001 and 
Francis, 2005) find a break in the absolute magnitude distribu- 
tion at Hj ~ 6.5. Indeed Francis (2005) finds that at higher 
absolute magnitudes the differential absolute magnitude distri- 
bution is virtually flat i.e. the cumulative increases linearly with 
Hj rather than as an exponential. From their analysis of LPCs 
that come close to Earth Fernandez & Sosa (2012) argue that 
the break in the absolute magnitude distribution is caused by a 
corresponding change in the slope of the size-frequency distri- 
bution rather than observational incompleteness. Thus in what 
follows we shall focus only on LPCs with Hj < 6.5. 

Everhart (1967) finds that 8 000 comets with H T < 10.9 
and q < 4 AU should pass through perihelion in 127 yr. He 
uses a photometric index v — 4 so that for his data Hj = H\q. 
From his paper the flux ratio between comets with Hj < 6.5 
and H T < 10.5 is 12, so that his flux for LPCs with H T < 6.5 
is only 3.3 yr~' for q < 2.5 AU. Approximately a third of 
LPCs are new (Wiegert & Tremaine, 1999; Fernandez & Sosa, 
2012), and thus Everhart's flux of new comets is 1.1 yr~' with 
q < 2.5 AU. Hughes (2001) quotes an LPC flux of 0.53 yr 1 per 
unit perihelion with Hj < 6.5 and he also assumes a photomet- 
ric index v — 4. When only a third of these are new his flux with 
q < 2.5 AU becomes 0.44 yr -1 . Francis (2005) determines the 
flux of new comets to be 2.9 yr~' with Hj < 10.9 and q < 4 AU. 
For new comets he uses a photometric index v = 2.44 on the 
inbound leg, which is Whipple's (1978) average for NCs. Ap- 
proximately 40% of his comets with Hj < 10.9 have Hj < 6.5 
and thus his flux of new comets is 0.73 yr -1 . Taking the average 
of the three sources, we arrive at an LPC flux of approximately 
(0.76 + 0.33) yr 1 with q < 2.5 AU and H T < 6.5. We shall 
use this to determine Noc below, but first we need to determine 
the size of an LPC with Hj < 6.5. This is done in the next 
subsection. 

5.2. The nuclear absolute magnitude of LPCs with Hj — 6.5 

There have been a number of attempts to relate the cometary 
absolute magnitude to the absolute magnitude of the nucleus, H 
(Bailey & Stagg, 1988; Weissman, 1996). However, these are 
generally unreliable. It is now well recognized that cometary 
nuclei develop non-volatile, lag-deposit crusts that reduce the 
fraction of the nucleus surface available for sublimation (Brin 
& Mendis, 1979; Fanale & Salvail, 1984). For most JFCs, 
the 'active fraction', that is the fraction of the nucleus surface 



area that must be active to explain the comet's water produc- 
tion rate, is typically only a few per cent, or even a fraction of 
a per cent (e.g. Fernandez et al., 1999). For LPCs, however, 
the active fraction is very large, and can 'exceed' 100% (Sosa 
& Fernandez, 2011). Thus LPCs are brighter than JFCs of a 
comparable size at the same heliocentric distance. Given these 
higher activity levels, the discovery probabilities for LPCs with 
small nucleus sizes should be considerably higher than those 
for JFCs of the same size. If the typical JFC is larger than 2 km 
in diameter (the median value for 67 measured JFC nuclei is 
3.7 km; Snodgrass et al., 201 1), it is entirely likely that an LPC 
with comparable brightness has a smaller, possibly sub-km nu- 
cleus. Is there some way of comparing values of Hj versus HI 

Sosa & Fernandez (201 1) use observational data of the wa- 
ter production of LPCs and the subsequent non-gravitational 
forces to derive a relation between the diameter of the nucleus 
of the LPC and its total absolute magnitude, which is given by 

logD = 1.2-0. 13H T , (4) 

where D is given in kilometres. Equation © appears valid 
mostly for small comets (Fernandez & Sosa, 2012). Equation 
(0 should be compared with the usual equation relating diam- 
eter and absolute magnitude for non-active objects, for which 
logD oc 0.2//. This difference in slope implies that for LPCs 
H - Hj is not a constant, but depends on H itself. 

Substituting Hj = 6.5 into equation © yields Dlpc = 
2.3 km, which corresponds to a nuclear magnitude of H = 17.3 
for a comet with a typical albedo of 4% (Fernandez et al., 2001). 
However, the above relation does not hold for JFCs, because 
these comets typically have a lower activity level than LPCs 
(e.g. Fernandez et al., 1999; Tancredi et al., 2006). To compare 
apples to apples we need to know what is the value of Hj for a 
JFC with H = 17.3. 



5.3. The total absolute magnitude of JFCs with H — 17.3 

Fernandez et al. (1999) plot the value of H - Hj as a func- 
tion of perihelion distance and fraction of active surface for 
JFCs with radii between 1 km and 5 km. Limiting ourselves 
to JFCs with q < 2.5 AU, we find H T ~ 9 when H ~ 17 
from their scatter plots, and thus the typical difference is 8, give 
or take 1 magnitude. Fernandez & Morbidelli (2006) also find 
H - H T - 8 ± 1 from small, faint JFCs with q < 1 .3 AU. In what 
follows, we consider a JFC with H = 17.3 to have Hj = 9.3, 
with an error of about 1 magnitudes. In conclusion, we find 
that LPCs with Hj < 6.5 have the same nuclear magnitude (i.e. 
physical size) as JFCs with Hj < 9.3, give or take a magni- 
tude. Thus, the OC to SD population ratio has to be determined 
from LPCs and JFCs with Hj < 6.5 and Hj < 9.3 respectively, 
rather than from populations with the same limiting total abso- 
lute magnitude. We do this below. 
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5.4. The Oort cloud to Scattered Disc population ratio 

To calibrate the number of SDOs with Hj < 9.3 and subse- 
quently compare this to the number of comets in the OC with 
the same diameter, we need to know the number of visible JFCs 
with Hj < 9.3. From their numerical simulations including 
physical ageing effects Di Sisto et al. (2009) find that there are 
117 active JFCs with q < 2.5 AU and D > 2 km. Levison & 
Duncan (1997) estimate the number of JFCs with q < 2.5 AU 
and Hj < 9 is 108, similar to that derived by Di Sisto et al. 

(2009) . 

The uncertainty in H - Hj for the JFCs is of the order 
of 1 magnitude. The cumulative absolute magnitude distri- 
bution of the comets obeys N(> H) oc \0~ ali . This corre- 
sponds to a cumulative size-frequency distribution of JFC nu- 
clei A^(> D) oc D y where y = 5a. For JFCs with diameters 
between approximately 2 km and 10 km the slope y ~ 2 (e.g. 
Lowry & Weissman, 2003; Meech et al., 2004; Snodgrass et al., 
201 1), corresponding to a — 0.4. With this value of a the JFC 
population varies by a factor of 2.5 for every magnitude differ- 
ence between H and Hj. Approximating the errors as Gaussian 
we have A^ v jfc = 1 17 + 50 and the median total number of bod- 
ies in the SD with D > 2.3 km is then N SD = 1.7^ x 10 9 . This 
is only a factor of three higher than Duncan & Levison (1997). 

We compute the number of comets in the OC with H < 17.3 
from the flux of new comets. As discussed earlier, this flux is 
about (0.76 ± 0.33) yr 1 for LPCs with H < 17.3 and q < 
2.5 AU. Kaib & Quinn (2009) state that the average fraction of 
the OC that has objects on orbits with q < 3 AU is 10~ n yr _1 . 
Thus the total OC population for comets with H < 17.3 is then 
Noc = (7.6 + 3.3) x 10 10 . Once again assuming the OC pop- 
ulation is normally distributed, the median population ratio be- 
tween the OC and SD for objects with D > 2.3 km is then 44+34. 
This is a factor of four higher than the results from our simula- 
tions presented in the previous section. However, the error bar 
overlaps with the nominal value from the simulations. 

What are the uncertainties in the above estimate? The above 
estimate of the OC/SD population ratio has taken into account 
the uncertainties in the SD decay rate, visible JFC production 
and H - Hj for JFCs and LPC flux. We did not yet consider 
other errors in the OC population. 

The most reliable estimates of the population of the OC 
yield of the order of 10 11 comets, but it is likely only to be cor- 
rect within a factor of two, or possibly even lower (Neslusan, 
2007). From the analysis of the motion of 26 LPCs and ac- 
counting for non-gravitational forces Krdlikowska & Dybczynski 

(2010) argue that approximately half of the LPCs that tradi- 
tionally were designated as being 'new' may in fact already be 
'old'. This result would yield a revised population estimate of 
the OC (4 x 10 10 ) that brings the OC to SD population ratio to 
23+jj, the error bar once again overlapping with the nominal 
value. 



In summary, we find that our scenario of the contemporary 
formation of the OC and SD in the framework of the Nice pre- 
dicts a OC/SD ratio that is about four times lower than the ra- 
tio deduced from observations. The latter, however, has large 
uncertainties, but the model-predicted and observation-deduced 
values do agree within the error bars. Therefore, the agreement 
between the simulations and observations cannot be rejected by 
a null hypothesis and thus there might be no problem with our 
scenario of the formation of these two comet reservoirs. 

6. Conclusions 

We have performed simulations of the formation and evolu- 
tion of the OC and SD in the framework of the Nice model. For 
OC formation simulations we kept the Sun in the current Galac- 
tic environment. The simulations lead to a somewhat higher 
capture efficiency than those of the more classical model where 
the giant planets are assumed to be on current orbits (e.g. Dones 
et al., 2004; Kaib & Quinn, 2008; Dybczynski et al., 2008; 
Brasser et al., 2010). We find that the efficiency of trapping 
comets in the OC is ~ 7% and that the simulated OC to SD 
population ratio is approximately 12 + 1, somewhat lower than 
most earlier results but still of the same order. 

We have shown that the SD produced in our simulations 
generates a population of JFCs that is consistent with the ob- 
served population. This is the first time that a dynamically hot 
SD is shown to be consistent with the JFC population. A pre- 
vious model of JFC origin started from a dynamically cold SD 
(Duncan & Levison, 1997), whose existence is challenged by 
observations. 

Using the link that we have established between the SD and 
the JFC population, as well as the link between the OC and the 
LPC population described in Wiegert & Tremaine (1999), we 
deduced a OC/SD population ratio from the observed fluxes of 
LPCs and JFCs. We performed the calculations for size-limited 
samples of comets (i.e. for LPCs and JFCs with the same di- 
ameter of the nucleus), by using the most recent conversions 
from total magnitude to nuclear magnitude available in the lit- 
erature for LPCs and JFCs. We found a population ratio of 
A^rat = 44+34, rou ghly a factor of four higher than the simulated 
nominal value of 12. The error bar of the observed ratio over- 
laps with the nominal value from simulations. This result takes 
into account all known uncertainties. Thus we conclude that 
our scenario of contemporary formation of the OC and SD in 
the framework of the Nice model is not inconsistent with the 
observations. 

Our scenario has several implications. First, given that on 
average the current SD population retained just 0.95% of the 
original trans-Neptunian disc population, we can estimate that 
the latter contained 1.9 x 10 11 comets with H < 17.3 (D > 
2.3 km) at the time of the instability of the giant planets. If 
instead we use the current OC population for this estimate, we 
find 10 12 comets. For comparison, the model by Morbidelli 
et al. (2009) of the primordial trans-Neptunian disc predicted 
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5x10 comets for the same value of H. 

Second, in our scenario both the OC and the SD are de- 
rived from the same parent population, i.e. the primordial trans- 
Neptunian disc. Thus, the LPCs and JFCs that come from these 
reservoirs should share (on average) the same physical proper- 
ties. This means that they should have the same size distribu- 
tion and the same range of chemical compositions. Both the 
size distributions and chemical compositions are very uncertain 
at the current stage of observational art. However, we remark 
that the size distribution of LPCs brighter than Hj = 6.5 and 
that of JFCs brighter than Hj — 9 are fairly similar; both are 
compatible with a cumulative size distribution with a slope of 
y = -2 (Lowry & Weissman, 2003; Meech et al., 2004; Snod- 
grass et al., 2011; Fernandez & Sosa, 2012). We also remark 
that the slope of the size distribution is very different for fainter 
JFCs, but this may be due to observational incompleteness or 
break-up of the comets on their way into the inner solar system. 
Regarding the chemical compositions, A'Hearn et al. (2012) 
argue that, from a statistical point of view, LPCs and JFCs are 
indistinguishable. Thus, although the last word has still to be 
said from the observational viewpoint, the prediction provided 
by our scenario seems to be verified. 
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